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Abstract 

f^yifferential cross sections of elastic neutron-proton scattering 
^3ia-ve been measured for the energy range from 200 MeV to 580 
The angular interval for the detection of the recoiling pro- 
ton ranges up to 48° in the laboratory system. This corresponds 
^Tpi an interval of the scattering angle from 80° to 180° in the 
^«entre-of-mass system. For absolute normalization the simulta- 
t~-aeously measured np — > d7r° reaction was used above 280 MeV. 
O^he charged ttNN coupling constant has been determined to 
l^np = 0.076 ± 0.001. 



Introduction 



^^eutron-proton elastic scattering at backward angles in the 
f^^edium energy regime has been the subject of several exper- 
i^Dments [1-6]. The common feature of these investigations is 
gteep rise of the differential cross section towards the back 
c^catterinp- angle of 180°. The slope of the sharp backward 
C^ak suggests a connection to the one pion exchange (OPE) 
f^m plitude. Several suggestions have been made for a the- 
(U)retical description of the experiments [7-16]. All of these 
i-^^roposals are able to describe the backward spike, but most 
them are purely phenomenological and some of them fail 
to describe other observables. 



^ .Precise backward scattering data offer the opportunity to 
C^valuate the pion-nucleon coupling constant /^nn [0 . The 
medium energy region is particularly suited for this pur- 
pose, because the pion pole is not very far from the physical 
region, so that one could expect a reliable extrapolation to 
it. However, the values determined so far by this method 
are not all in accordance with each other and with the value 
obtained from pion-nucleon scattering [Q . 

The results presented in this paper have been obtained 
from four separate experiments of different angular ranges, 
labeled I-IV [4, 19-21]. Together they span the interval from 
about 80° to 180° in the centre-of-mass system. The neu- 
tron energies range from 200 MeV to 580 MeV in steps of 
20 MeV. 



2 Experiment 

The experimental set-ups and techniques of experiments I- 
IV have been quite similar but with differences in detail of 



the accelerator performance, the beam arrangement and the 
detection equipment. Here we outline the common features 
only and refer to a forthcoming paper 1 19 for details. 

The experiments have been performed at the Paul Scher- 
rer Institute (PSI), the former Swiss Institute for Nuclear 
Research (SIN). The proton beam of the ring cyclotron of 
589 MeV energy consists of bunches with a width of less 
than 1 ns at a rate of 50.63 MHz for experiments I and 
III and 16.88 MHz for II and IV, corresponding to bunch 
spacings of 19.75 ns and 59.25 ns, respectively. The beam 
current during the data taking was 60-100 /xA. The neu- 
trons were produced on a thick target of beryllium (I-IH) 
or carbon (IV). Neutrons escaped through a collimator hole 
in the beam dump at an angle of 60 mrad with respect to 
the incident protons. They were shaped by two additional 
collimators to a beam of about 2x2 cm^ at a distance of 
61m from the production target. Charged particle contam- 
inations were eliminated by cleaning magnets behind the 
collimators. A lead filter was inserted in order to reduce 
the photon component of the beam which mainly originates 
from the decay of neutral pions in the production target. 

The continuous neutron energy spectrum consists of the 
40 MeV wide quasielastic peak at about 540 MeV and a 
broad distribution at lower energies, which can be ascribed 
to inelastic processes [^. The detailed shape of the spec- 
trum depends on the target material and its thickness. 

After a flight path of 61 meters, kept at rough vacuum 
of about 100 Pa, the neutron beam hits a liquid hydrogen 
target. For the analysis of the scattering products a magnet 
spectrometer was installed on a turn table. The spectrome- 
ter was equipped with drift chambers and two scintillation 
counters, a thin one (1 mm thick) in front of the magnet, 
and a thicker one (1 cm thick) behind it. They allow to 
measure both, the time-of-flight of the particle detected by 
the spectrometer, and that of the incoming neutron with 
respect to the rf-signal of the accelerator. The angular ac- 
ceptance is almost 20°. The average momentum resolution 
is about 3 % FWHM. 



3 Data taking and analysis 

The data have been taken at different run periods, each 
extending over several weeks. The contribution of the tar- 
get surroundings and spectrometer materials was measured 
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with an empty target. It was subtracted after normaliza- 
tion to the neutron intensity of the full target measurement. 
Corrections have been applied for the event rate dependent 
dead-time losses. 

The long periods of data taking required special attention 
to the long term variations of the whole system, particularly 
of the neutron time-of-flight measurement. The stability of 
the electronics has been checked regularly and close control 
of the drift chamber gas flow and high voltage supply has 
been maintained. 

The neutron intensity was monitored in three different 
ways. The integrated primary proton beam intensity was 
provided from the accelerator control centre. This signal, 
though convenient, was not sufficiently reliable because of 
variations of the focusing of the primary beam on the pro- 
duction target. The second monitor consisted of a three 
stage scintillation counter telescope, which recorded charged 
particles emerging from a thin polyethylene target placed in 
the beam at about 44 m from the neutron target. A third 
monitor was installed behind the magnet spectrometer. It 
recorded elastic np scattering events from a polyethylene 
block in coincidence. 

These monitors were intended to measure the relative 
neutron intensity. Above 280 MeV the absolute normaliza- 
tion was performed by the simultaneously recorded np — > 
d7r° reaction as reference cross section. This is discussed in 
sect. p. 

3.1 Time-of-flight calibration 

As mentioned above, the energy spectrum of the incident 
neutrons is continuous. Therefore, the neutron time-of- 
flight measurement is the basis of the incident energy de- 
termination, and the control of its stability is crucial. It 
requires careful calibration of the zero point and the conver- 
sion gain of the time-to-digital converters. For the control 
of the time zero point we used the high energy photons in 
the beam. The lead filter in the neutron beam was removed 
for these calibration runs, and the hydrogen target was re- 
placed by a lead slab for a higher conversion rate. With 
the spectrometer field inverted and reduced appropriately, 
the converted electrons were detected. From the width of 
the peak the time resolution has been evaluated to be 0.8 
ns FWHM including the bunch width of the primary beam 
and the contribution from electronics. The corresponding 
neutron energy resolutions vary from 1.2 MeV at 200 MeV 
to 7.6 MeV at 580 MeV. 

These calibrations have been performed in regular inter- 
vals. Deviations from the overall mean have been corrected 
for each run. 

3.2 Data reduction and event selection 

The first step in the off-line analysis was a reduction of the 
data by setting cuts in order to select good events. These 
cuts included a unique track in the drift chambers, start- 
ing in the target volume and emitted within the full accep- 



tance of the spectrometer in vertical and horizontal direc- 
tion. Next, the mass of the particle in the spectrometer 
was determined from the measured time-of-fiight through 
the spectrometer and the measured momentum. A mass 
selection has been applied by rejecting masses m < 0.5 nip 
and m > 1.5 mp, where mp is the proton mass. The up- 
per limit was increased to m > 2.5 mp for experiments I 
and II, in order to keep also deuteron events, which were 
used for the normalization of the np data (see sect. P). Fi- 
nally, low momentum particles have been excluded by cuts 
at p = 500 MeV/c for protons and p = 620 MeV/c for 
deuterons. 

3.3 Neutron energy determination 

The emission angle of the recoil proton is determined by 
the set of drift chambers in front of the magnet, and its 
momentum by the deflection in the magnetic field. The 
proton energy is obtained after correction of the energy loss 
in traversing the material of target and spectrometer. The 
incident neutron energy is determined from the measured 
neutron time-of-flight ti. A complication arises from the 
fact, that the bunch interval given by the rf-signal (19.75 
ns or 59.25 ns) is much shorter than the time-of-flight from 
the production target to the hydrogen target (more than 
200 ns) . This introduces ambiguities on the proper multiple 
of the bunch spacing time which has to be added to the 
measured value. However, by comparison with the time- 
of-flight icaic,£, calculated from the measured recoil proton 
momentum vector, the correct number of bunch spacings 
can be determined. 

4 Experimental results 

4.1 np difl"erential cross sections 

The data were binned in energy intervals of 20 MeV, and 
angular bins of 0.5 degrees in the centre-of-mass (CM) sys- 
tem. Each bin has been corrected for absorption of the 
recoil protons in the target and the spectrometer material. 
This energy dependent correction was below 1 % in all cases. 
Another correction was necessary due to a small polariza- 
tion component of the neutron beam in vertical direction 
of .05 < |Pn| < -08. This comes about since the neutrons 
are produced at an angle of 60 mrad (cf. sect. |^). The en- 
ergy and angle dependent effect of the analyzing power of 
the np elastic scattering has been taken into account. The 
correction to the differential cross section was below 2.5%. 
The errors of these corrections contribute to the errors of 
the differential cross sections with less than 0.3%. 

The four sets of data have been combined in an appro- 
priate way. For the experiments I and II a common angular 
region exists which has been used to evaluate the multiply- 
ing factor for experiment I. The uncertainty introduced by 
this procedure for the different energies has been determined 
as 0.6% from the spread of the individual factors from the 
average. In a similar way the data set from experiment III 



2 



200 MeV ^ 
i 

*»«»>«****** 


220 MeV / 


: 240 MeV / 

/ 


: 260 MeV / 

^ .y 

********* 
**«*.********^ 


: 280 MeV ^ 
i 

■ **************** 


: 300 MeV / 
* 
/ 

»♦** 

. ««.*«»*«*"**^ 


: 320 MeV / 
/ 

y 


: 340 MeV / 

* 
/ 

****** 
. *****************^ 


360 MeV / 

/ 


380 MeV ^ 

# 

. ********************* 


" 400 MeV Ji 
i 
/ 

y 


" 420 MeV ^ 
t 

/ 

**********************'** 


440 MeV / 

/ 

^ y 


460 MeV / 

* 
• 

/ 

y 


: 480 MeV / 

• 

* 
< 

^y^ 


: 500 MeV / 
t 

9 
/ 

y 


520 MeV i 

8 

J 


540 MeV 1 

» 

J 


: 560 MeV / 
# 
/ 


: 580 MeV ,* 

n^y'^ 



90 120 150 180 90 120 150 180 90 120 150 180 90 120 150 180 



Figure 1: Relative differential np cross sections as a function of the CM-scattering angle ^cm- °: data set I, 150° < ^cm < 180°; 
• : data set II, 140° < 6>cm < 180°; x: data set III, 129° < Qcw < 156°; *: data set IV, 80° < ^cm < 126°. The data are normalized 
to 1 at 6cM = 179.2°. Errors are not visible within the symbol size. 



has been linked to the combined data of set I and set II. 

Here the uncertainty was 0.7%. This procedure could not 
be followed for the normalization of experiment IV because 
no overlap of data points exists. We therefore fitted for each 
neutron energy the data points of experiment III for CM- 
angles smaller than 145° together with the same number of 
data points of set IV by Legendre polynomials of 4th order 
with a free scaling factor. Each data point has been given 
the same statistical weight of 2 % in the fit procedure in 
order to assign the same weight for the two data sets. The 
average value of for the relative normalization was 1.41, 
and the mean error of the scaling factor is 2 %. 

The statistical errors of the difi'erential cross sections 



range from 1 % to 2 %. The spread of the data with respect 

to smooth fits with polynomials of an appropriate order is 
more like 2-2.5 %. This reflects systematic uncertainties due 
to 

- short term changes in the phase of the rf-signal 

- imperfections of the drift chambers and other detectors 

- remaining effects of gain shifts of detectors and elec- 
tronics. 

Taken together they are in the order of 1-2 %, depending on 
the set-up. Systematic errors of 1.4%, 0.6% and 1.7% have 
been added in quadrature to the data sets II, III and IV, 
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respectively. A systematic error of 0.5 % had been added 
earlier to data set I p3[ . 

The relative differential cross sections for the four sep- 
arate experiments are displayed together in Fig. 0. The 
angular distributions are normalized to 1.0 at the largest 
angle of data set II. The numerical values are tabulated in 
ref. H). 

4.2 Discussion 

A common feature of the angular dependence at all ener- 
gies is the sharp rise towards the backward direction. It 
is followed by a less steep decrease passing through a wide 
minimum, which is shifted to smaller angles with decreasing 
energy. In the transition region of the two different slopes 
there is an indication of a bump. A similar shape is indi- 
cated in the LAMPF data j^] but not as clear. 

The energy dependence of the cross section at 180° is 
shown in Fig. ^. It is roughly constant in our energy range, 
and is well described by the phase shift predictions. With 
decreasing angle a deviation of the data from the phase shift 
predictions develops within the full energy range. This can 
be clearly seen by a comparison of the cross section ratios 
cr(6')/cr(180°) as shown in Fig. | for 9 = 135° and 90°. The 
deviations occur for both phase shift solutions of Arndt et 
al. ||2^ and Bystricky et al. ||2^. It reflects the fact that 
for the extreme backward scattering angles several precise 
and consistent measurements exist which pin down the 
phase shift solutions, whereas the angular region around 
= 90° has been covered only scarcely at singular ener- 
gies so far. The cross section ratios, as displayed in Fig. ^, 
are decreasing with increasing energy, except for the ratio 
cr(90°)/(T(180°) which is almost energy independent above 
250 MeV. 

5 Absolute normalization 

For the absolute normalization of the cross section the in- 
cident neutron intensity has to be known as a function of 
energy. In our case the neutron intensity was obtained for 
experiment II from a comparison with the simultaneously 
measured reaction 



np 



d^" 



(1) 



where the deuterons have been recorded by the spectrometer 
like the recoiling protons of the elastic scattering. 

Isospin independence is used to relate the cross section of 
the process (0) to the cross section of reaction 



PP 



dTT^ 



(2) 



for which precise data exist. However, isospin symmetry, 
which relates the two reactions by adTr° — ^0"^^+ is not 
exact since the masses in the initial and in the final states 
are different for (|l|) and (||). Also, the Coulomb interaction 
for process (||) must be taken into account. This is discussed 
in detail in section b3. 
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Figure 2: Comparison of our data with phase shift predic- 
tions. Full lines: Arndt et al. dotted lines: Bystricky et 
al. |2^]. In the upper diagram the cross sections extrapolated 
to 180° are shown. In the lower diagram the cross section ratios 
^(61)/ ^(180°) vs. Tn are shown for 6 = 135° and e = 90°. 



5.1 Measurement of ^(np — > d7r°) 

The angular acceptance of the magnet spectrometer is sufh- 
ciently large to cover the full angular range of the deuterons, 
^max ^ 22°, so that a complete angular distribution of reac- 
tion (d) is obtained. The energy of the deuterons emitted in 
([l]) varies strongly with the emission angle. For an incident 
neutron energy of 400 MeV, as an example, the deuteron 
energy varies from 130 MeV to 260 MeV. In this energy re- 
gion the deuteron cross section on nuclei is strongly energy 
dependent and distorts the angular distribution. Besides 
the energy loss the absorption by the spectrometer material 
has to be corrected for. While the energy loss correction is 
straightforward, the correction for absorption loss is more 
involved because of the missing knowledge of deuteron cross 
sections on nuclei. 

For the deuteron total cross section on a nucleus with 
mass number A we use [ETI 



adA = 0.97(o-nA + CFpA) ■ 



(3) 



The material around the target and the spectrometer con- 
sists mainly of hydrogen and the self-conjugate nuclei C, N, 
O and Ar, for which we assume a^A — CpA- The deuteron 
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cross section eq. (g) reduces then to 

(TdA = l-94(TnA 



(4) 



Neutron total cross sections on nuclei have been measured 
in a previous work |^^ . For N and Ar, which have not 
been measured directly, the cross section has been evaluated 
from a parameterization of the energy and mass number 
dependence given in . For the deuteron cross section on 
hydrogen adp at the energy T we used Cpd data or the (T„p 
data of H at T/2. 

The absorption correction is then obtained by summing 
up the contributions of the different elements according to 
their thickness in the spectrometer. The resulting correc- 
tion varies from 5.2% at Td = 120 MeV to 2.5% at 400 
MeV. As in the case of elastic np scattering (cf. sect. 4.1) 



the effect of a small polarization of the incoming neutrons 
has been corrected for. This angular and energy dependent 
correction for the deuteron production is below 0.6% for 
Tn < 490 MeV. At our highest energy it reaches 2.6%. The 
errors of both corrections taken together contribute to the 
errors of the differential production cross sections with less 
than 0.5%. 

5.2 Mass- and Coulomb corrections 

Niskanen and Vestama studied the symmetry violat- 
ing effects for the reactions (|l|) and (||) in the framework of 
a coupled-channel method with different pion production 
mechanisms. Because the final momentum is most rele- 
vant in threshold reactions with a large negative Q-value 
the comparisons are made at the same final pion momenta 
Ptt in the CM system, commonly expressed by the dimen- 
sionless variable 77 = pTr/mT^+. For both reactions, the pion 
production cross section can be factorized as dcr/dO = P R 
with the phase space factor P and R the sum of the squared 
matrix elements of the pion production operator. With this 
ansatz the relative change in the cross sections was deter- 
mined in good approximation as 



5(j _ 5P 
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(5) 



with 



da = 2 CTdTr" - Crd^+ (Tip 

SR = RdTT" — RdTT+ Ra- 
sp = Pd.o - Pdx+ Pa- 



i • (2o'd7r° + 0'd7r+) 
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Niskanen and Vestama calculated the effect due to the mass 
differences and incorporated the Coulomb correction in the 
final state only. At higher proton energies, the Coulomb 
correction for the pp initial state is comparable with the 
one in the d7r+ final state. Therefore we use the results 
of Rcf. for the change in the matrix element without 
Coulomb correction and apply Coulomb penetration factors 
for a point proton with an extended proton and for a point 
pion with an extended deuteron. The weighted sum of the 



Coulomb penetration factors for the different orbital angular 
momentum states with weighting factors from is used 
as Coulomb correction C^^ for the squared matrix elements 
i?Q(j^+ without Coulomb correction, 



R. 



'd7r+ 



•0d7r+ ■ 



(6) 



Table || shows some relevant information used to transform 
the cross section crd7r+ of (||) to (Td7r= of the process (|]) by 
the relation 



_ 1 



2 + Sa/a^ 
2 — Sa/ffa 



(7) 



The uncertainty for SR/R^v is according to Ref. |Q about 
0.01. Together with the Coulomb corrections which we ap- 
plied to the pp initial state, this increases to 0.015. This 
error contribution is included in the errors given in col. 9 of 
Table |. 

5.3 Parameterization of cr(pp — >■ d7r+) 

Cross sections for reaction (||) as well as for the inverse reac- 
tion TT+d pp have been measured with sufficient accuracy 
in the relevant energy range from threshold up to 640 MeV. 
The data base [31-35] consists of 61 integrated cross sections 
for the reaction between 288 and 641 MeV incident pro- 
ton energy, and of 74 cross sections for the inverse process 
between 1.8 and 174 MeV incident pion energy. The latter 
ones have been converted to (j^) via detailed balance. Sta- 
tistical and systematic errors of each data point have been 
added in quadrature (Gaussian) . The data set has been fit- 
ted by a parameterization similar to Bystricky et al. 
as a function of rj. The fit error of adTT+ ranges from 2.1 % 
at ?7 = 0.396 to 0.7% at r? = 1.582. 

5.4 Absolute normalization above 300 MeV 

The absolute normalization factor for the elastic np scatter- 
ing data is given with respect to the highest measured angle 
6»„iax = 179.2° of data set II 



da- , , _ N^piOmax) 

diZ j Nd^odi2d7r° 



(8) 



where A^np and N^^o are the solid angle corrected event 
numbers for the two processes in corresponding energy bins. 
The results are given in col. 10 of Table |. Close to pion 
threshold and below, this method can not be applied. 



6 The ttNN coupling constant / 
6.1 The Chew method 
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ttNN 



The value of the ttNN coupling constant /^nn focussed 
new interest recently in connection with a discussion on a 
conceivable breaking of charge independence [37-41]. The 
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300.2 
320.1 
340.0 
360.1 
379.9 
400.2 
420.3 
440.5 
460.4 
480.5 
501.0 
520.8 
540.3 
559.6 
578.1 



0.3960 
0.5369 
0.6527 
0.7558 
0.8486 
0.9375 
1.0206 
1.1004 
1.1760 
1.2498 
1.3229 
1.3916 
1.4577 
1.5217 
1.5820 



312.0 
331.4 
350.9 
370.7 
390.1 
410.1 
429.9 
449.9 
469.6 
489.5 
509.8 
529.4 
548.7 
567.9 
586.2 



0.9364 
0.9471 
0.9525 
0.9559 
0.9584 
0.9606 
0.9622 
0.9638 
0.9651 
0.9665 
0.9678 
0.9690 
0.9700 
0.9711 
0.9719 



0.1194 
0.0950 
0.0778 
0.0643 
0.0521 
0.0412 
0.0325 
0.0254 
0.0200 
0.0146 
0.0132 
0.0225 
0.0400 
0.0560 
0.0700 



-0.0095 
-0.0084 
-0.0074 
-0.0066 
-0.0059 
-0.0053 
-0.0048 
-0.0043 
-0.0039 
-0.0036 
-0.0032 
-0.0030 
-0.0027 
-0.0025 
-0.0023 



0.1099 
0.0866 
0.0704 
0.0577 
0.0462 
0.0359 
0.0277 
0.0211 
0.0161 
0.0110 
0.0100 
0.0195 
0.0373 
0.0535 
0.0677 



0.1225 
0.2288 
0.3639 
0.5296 
0.7213 
0.9474 
1.200 
1.481 
1.777 
2.085 
2.392 
2.664 
2.889 
3.054 
3.143 



0.0684 
0.1248 
0.1952 
0.2805 
0.3777 
0.4910 
0.6171 
0.7563 
0.9030 
1.054 
1.208 
1.358 
1.499 
1.611 
1.682 



Table I: Absolute normalization of the elastic np cross sections. All cross section related quantities are given in the CM system. 
Col. 1: incident neutron energy; col. 2: 77 = p.^/m^+ for reaction col. 3: incident proton energy for reaction at same rj; col. 4-7: 
Contributions to Mass and Coulomb correction (details see text); col. 8: integrated cross section of reaction S); col. 9: integrated 
cross section of reaction dlh; col. 10: absolute normalization factor for np elastic cross sections 



charged coupling constant, can be obtained by an ex- 6.2 Test with pseudo data 



trapolation of the measured backward elastic differential 
cross section to the pion pole. The standard procedure in 
the past has been the Chew method |l^, 42 1, which ex- 
trapolates the Chew funcion y{x) defined by 



y('^) (-^ -^-K—pole)" 



df^ 



J=0 



(9) 



to the pion pole in the unphysical region. Here x — cost 
with 9 the CM-angle of the scattered neutron and 



•^Tz—pole 



COS 9 



TT—pole 



2E,,Ep)/2pl (10) 



with 



and nip the masses of the charged pion, the 



neutron and the proton, En and Ep the CM-energies of the 
neutron and the proton and Pn the CM-momentum of the 
neutron. The Pj are polynomials generated with special 
recurrence relations, so that they are orthogonal in the range 
where data points exist [Q. At the pion pole the Chew 
function gives 



y {•^77— pole) 

where s is the total energy squared. 



(11) 



We have tested this extraction method with pseudo data 
between 240 MeV and 540 MeV neutron laboratory kinetic 
energy in 100 MeV steps. The data were generated from 
the regularized OPE model of Gibbs and Loiseau jl^ in 
the angular range from 80° to 180° in 1° steps. The input 
coupUng constant was fixed at = 0.076. Uncertainties 
of the pseudo data were generated with a Gaussian random 



















0.08 




A| 


A 


° + 


6 








0.07 






+ 










0.06 
















0.05 


1 


1 






1 


1 


1 




4 


5 


6 


7 


8 


9 


10 



polynomial order 

Figure 3: fits at Tn — 440MeV vs. maximum polynomial 
order. Circles represent fits without, triangles with conformal 
mapping. Filled symbols are used for pseudo data with realistic 
errors (2%), whereas open symbols are used for high precision 
pseudo data (0.1%). 

error distribution. If these uncertainties are < 0.1%, the 
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model coupling constant can be reproduced with an error 
of « 0.25%. For pseudo data with uncertainties of 2%, 
corresponding to the present experiment, the error for 
is 2.4%. Moreover, the extracted coupling constant is sys- 
tematically smaller than its input value. This is because 
the extracted rises with the maximum polynomial order 
of the fit as shown in Fig. ^ (open circles). Applying an 
F test [Q to the fit gives a lower maximum order of the 
polynomials of eq. ^ for data with higher uncertainties. At 
Tn = 440 MeV, for instance, and for uncertainties < 0.1% 
the F test results in the maximum order jmax ~ 10 and 
/2 = 0.0761 ± 0.2%. For uncertainties of 2% the F test 
results in j^ax = 7 and = 0.0746 ± 2.4 %. 

6.3 Chew method with conformal mapping 

The convergence properties of the polynomial expansion is 
improved by conformal mapping of the variable x before ex- 
trapolating to the pole ji^. Wc used this method similar 
to Dumbrajs et al. |4^. The position of the neutral pion 
pole and the onset of the two pion (7r°7r+) exchange cut are 
mapped to (-1-1,0) and (-1,0) in the complex plane, respec- 
tively. The other exchange cuts lie on the unit circle. In 
addition, the data points are symmetrized around zero. 

This mapping method was tested as before with pseudo 
data. Again, the model coupling constant is reproduced 
for uncertainties < 0.1%. The convergence is faster than 
without conformal mapping, and the extrapolation error is 
of the same order or slightly larger. For the example at 
Tn = 440 MeV given above, the F test results a value of 
/c. map — 0.0759 ± 0.3% with a reduced polynomial order 
of jmax = 8 for uncertainties of < 0.1%. With uncertain- 
ties of 2 %, the F test results in the lower value jmax = 5 
and /c map — 0.0786 ± 1.6%. Thus, the conformal map- 
ping method results in a systematic upward shift for the 
extracted coupling constant, as shown by the triangles in 
Fig. ^. Again, this comes about because the extracted 
/c map depends on the maximum polynomial order, which 
is lower for less precise pseudo data. 

Similar results were obtained with pseudo data from phase 
shift predictions of Arndt et al. [p5| . Again, the input 
value = 0.076 could be reproduced only with high pre- 
cision (< 0.1%) pseudo data by both methods. For uncer- 
tainties of 2 % the result was systematically too low with- 
out, and too high with conformal mapping. Thus, both 
methods have systematic offsets of opposite sign. In prin- 
ciple, they can be taken into account by a correction factor 

^corr = /hiput/ /extracted- 

6.4 Determination of f ^ from our data 

Besides the systematic uncertainties caused by the errors of 
the data, there is also a dependence on the angular range 
used in the fit. With decreasing range also the extracted 
coupling constant is decreasing. This can be seen from the 
results on and /^^map obtained with the real data at 
Tn = 440 MeV, as shown in Fig. ||. For the fits at this 
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Figure 4: Dependence of fits on the angular range at Tn = 
440 MeV. Full dots are the results with conformal mapping. 

energy all data points were given a statistical error of 2.8 %. 
With this choice the reduced of the fit is of the order 
of 1 for the whole angular range. The error of is about 
15% for an angular range A6 = 15° and drops to about 
3% for A6 = 100°. According to Fig. ^ an angular range 
A9 > 50° is needed. The variation of the extracted 
values is minimal for 75° < A9 < 95°. For all energies it was 
found, that the angular range should extend to the region 
where the differential cross section reaches its minimum. In 
order to have no different weights from the four data sets, 
we have chosen for each laboratory kinetic energy a common 
relative error per data point, which was adjusted to get a 
reduced Xr of the order of 1 for the whole angular range. 
The maximum polynomial order jmax was determined either 
when a minimum of was reached or when < 4, which 
corresponds to a significance of at least 95 % for the term 
with jmax- The fit error of is then below 5 %. The results 
of the fits at energies where an absolute normalization exists 
are collected in Table |l[ 

The correction factors /ccorr , which were needed to correct 
for the systematic offset introduced by the Chew extrapo- 
lation were determined for both methods by an extended 
study with OPE pseudo data. For each energy we have gen- 
erated 1000 sets of pseudo data, with randomly distributed 
errors s as given by the experimental data (see Table ^ . For 
the Chew method without conformal mapping no unique 
jmax was found for a given energy. Therefore we have used 
the fccorr obtained for that value of jmax which was obtained 
for the real data. This complication does not arise with 
conformal mapping, where the fit criteria (F-test) lead to 
the same value of jmax at all energies. Moreover, the energy 
variation of fccon- is only « 1 % as compared with almost 
10 % without conformal mapping. Therefore we prefer to 
extract an energy independent value of using the con- 
formal mapping method. The maximum polynomial order 
jmax, the resulting value of the fit error A/^ and the 
correction factor /ccorr is given in Table |l| for both methods 
at each energy. 
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input data 


results of Chev 
without conformal mapping 


f extrapolation 

with conformal mapping 


T„ 
(MeV) 


A9 

n 




s 
(%) 


Jmax 


/c 


A/2 

(%) 




Jmax 


J c, map 


A/2 

(%) 




300.2 


94.2 


108 


2.3 


7 


0.0766 


4.7 


0.980 


5 


0.0730 


3.1 


0.977 


320.1 


92.5 


106 


2.3 


7 


0.0778 


4.2 


0.983 


5 


0.0762 


2.8 


0.975 


340.0 


88.8 


106 


2.4 


6 


0.0713 


3.4 


1.040 


5 


0.0765 


3.1 


0.973 


360.1 


87.0 


106 


2.4 


7 


0.0794 


4.4 


0.985 


5 


0.0807 


3.2 


0.974 


379.9 


85.2 


104 


2.5 


7 


0.0791 


4.5 


0.982 


5 


0.0778 


3.2 


0.971 


400.2 


83.5 


104 


2.8 


6 


0.0713 


3.9 


1.049 


5 


0.0786 


3.6 


0.969 


420.3 


83.7 


106 


2.8 


6 


0.0700 


3.8 


1.056 


5 


0.0776 


3.4 


0.970 


440.5 


79.9 


103 


2.8 


6 


0.0742 


3.5 


1.049 


5 


0.0797 


3.6 


0.967 


460.4 


80.1 


104 


2.5 


6 


0.0725 


3.1 


1.060 


5 


0.0780 


3.1 


0.967 


480.5 


78.3 


103 


2.6 


6 


0.0734 


3.2 


1.060 


5 


0.0792 


3.2 


0.965 


501.0 


74.4 


101 


2.5 


6 


0.0757 


3.0 


1.053 


5 


0.0783 


3.3 


0.965 


520.8 


74.6 


104 


2.1 


6 


0.0769 


2.3 


1.060 


5 


0.0811 


2.6 


0.964 


540.3 


74.9 


104 


2.2 


6 


0.0746 


2.4 


1.063 


5 


0.0794 


2.6 


0.963 


559.6 


75.1 


103 


2.5 


6 


0.0741 


2.6 


1.077 


5 


0.0815 


2.7 


0.963 



Table II: Fit results for obtained with the standard Chew method (/c), and with conformal mapping (/c, map) for all energies 
with absolute normalization. A^: angular range; A'': number of data points; s: error of individual data point; jmax: maximum 
polynomial order according to F test; A/^: fit error of and /c, map! respectively; fecorr: correction factor for systematic offset. 



A contribution to the error which has not been taken 
into account so far is due to the uncertainty of ~ 2 % in 
the relative normalization of data set IV with respect to 
the other sets. It was investigated systematically and a 
variation of about 1.5% of the extracted /2 was observed. 
This contribution is not contained in A/2 of Table ||. It 
was, however, added in quadrature to each A/2 before the 
energy averaged weighted mean (/2„p) = {p ■ fccorr) has 
been calculated. 

For the prefered method with conformal mapping we ob- 
tain as our final result the weighted mean of the charged 
coupling constant in the energy range 300 MeV < Tn < 
560 MeV: 

Lflnp) = 0-0760 ± 0.0008 . 

The given error is the propagated error calculated from the 
individual A/2. It is slightly larger than the standard de- 
viation of the mean. 

It has to be kept in mind, that both, the relative error 
of the data points as well as the absolute normalization en- 
ter and limit the determination of /2. In our case both 
error types are of the same order. With more precise data, 
i.e. smaller errors of the relative cross sections, the fitting 
procedure would be more stable and the systematic effects 
be smaller. Changes in the absolute normalization factor 
and its error, on the other hand, propagate to the /2 values 
only with a factor 1/2 since the Chew function at the pole 
depends on /''. 

All method tests and the calculation of fccorr have been 
performed with OPE model data. Therefore, a model con- 
tribution should eventually be added to the error margin. 

Our value of 0.0760 for the charged coupling constant 



is remarkably lower than the value obtained by Hohler et 
al. jl^ from pion-nucleon scattering. It is, however, in good 
agreement with more recent determinations of /2 from both, 
pion-nucleon partial wave analysis [ ^7| and nucleon-nucleon 
scattering pTj . 
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